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In this article we address the problem of getting the temperature dependence of the tt-tt scattering 
lengths in the frame of the linear sigma model. Using the real time formalism, we calculate all the 
relevant one loop diagrams. The temperature corrections are only considered in the pion sector, 
due to the Boltzmann suppression for heavier fields like the sigma meson. From this analysis we 
obtain the thermal behavior of the s-waves scattering lengths ag and associated to isospin 7 = 
and 1 — 2, respectively. If we normalize with the zero temperature value it turns out that 



grows with temperature, whereas the opposite occurs with ° y . Finally we compare our results 

with other determinations of the scattering lengths based on the Nambu-Jona-Lasinio model and 
chiral perturbation theory. 
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I. INTRODUCTION 



The discussion of the thermal evolution of tt-tt scat- 
tering lengths turns out to be a relevant problem in 
the context of heavy ion collisions. In fact, we know 
that in the central rapidity region, precisely where the 
quark-gluon plasma is expected to be created, a big 
amount of thermalized pions will be produced. Those 
pions will interact among themselves and the tt-tt scat- 
tering lengths are crucial parameters in order to describe 
the scattering amplitudes. At zero temperature the 
TT-TT scattering lengths were first measured by Rosellet 
et al A recent review about the present status 

of the experimental measurements of tt-tt scattering 
lengths and their comparison with different theoretical 
approaches can be found in 2]. 

In this article we will present a detailed calculation of the 
thermal corrections to the tt-tt scattering lengths in the 
frame of the linear sigma model il. As it is well known, 
the linear sigma model is an effective, renormalizablc 
low energy description of hadron dynamics. Our 
calculations will be done using the real time formalism 
at the one loop level. 



II. LINEAR SIGMA MODEL AND tt-tt 
SCATTERING 



scalar field cr. The idea is to define a new field s such 
that a = s + V. Obviously (s) = 0. ip corresponds 
to an isospin doublet associated to the nucleons, tt 
denotes the pion isotriplet field and ca is the term that 
breaks explicitly the SU{2) x SU{2) chiral symmetry, 
e is a small dimensionless parameter. It is interesting 
to remark that all fields in the model have masses 
determined by w. In fact, the following relations are 
vaKd: = gv, = + '^^^^ and = ^'^ + 3A^w^. 
Perturbation theory at the tree level allows us to identify 
the pion decay constants as = v. This model has 
been considered in the context of finite temperature 
by several authors, discussing the thermal evolution of 
masses, U{T), the effective potential, etc Q, [1], 0, [1]. 



C = i'lij^df, - - g{s + in ■ T-/5)]ip + -^[{dnf + w^tt^] 
+ ^[{daf + mis'] - WKsis' + TT^) - ^(s'- + n'f 



+ {ec — vm1)s 



(1) 



Since our idea is to use the linear sigma model for 
calculating tt-tt scattering lengths, let us remind briefly 
the formalism. The scattering amplitude has the general 
form 



The linear sigma model in the phase where the chiral 
symmetry is broken is given by the lagrangian below, 
where v ~ (a) is the vacuum expectation value of the 
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Tafs-s-y = t, u)daf3Ss^ + A{t, S, u)5a^5i3S 

+ A{u,t,s)5as5l3^. (2) 

where the a, (3, 7, 5 denote isospin components. 

By using appropriate projection operators, it is pos- 
sible to find the following isospin dependent scattering 
amplitudes 
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r" = 3A(s, t, u) + Ait, s, u) + A{u, t, s) (3) 

T'^ = A{t,s,u) - A{u,t,s) (4) 

= A(t,s,M) + A(M,t,s), (5) 

where denotes a scattering amplitude in a given 
isospin channel. 

As it is well known Q , the isospin dependent scattering 
amplitude can be expanded in partial waves T/. 



1 

T/(s) = — y d{cos0)Pi{cos9)T'{s,t,u) (6) 



Below the inelastic threshold the partial scattering am- 
plitudes can be parametrized as [lO| . 



s — 4m7r^ 2i 



1), 



(7) 



where Sf is a phase-shift in the £ channel. The scattering 
lengths are important parameters in order to describe low 
energy interactions. In fact, our last expression can be 
expanded according to 



Experimental 
Results 



0.26 ±0.05 



0.25 ± 0.03 



-0.028 ±0.012 



-0.082 ± 0.008 



0.038 ± 0.002 



Chiral Perturba- 
tion Theory 



0.16 



0.18 



167r/j 



-0.044 



-0.089 



0.030 







Linear 
Model 



Sigma 



: 0.22 



1287r/. 



2- = 0.27 



-0.044 



-0.089 



4^ 



= 0.030 



48T7f 



: 0.015 



TABLE I: Comparison between the experimental values [l|, 
first order prediction from chiral perturbation theory [ll|] and 
our results at the tree level. 
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4, ,2 



t — mi 



t — mi 



4AS' 



4„2 



(10) 

(11) 

(12) 



(8) 




The parameters and are the scattering lengths 
and scattering slopes, respectively. In general, the scat- 
tering lengths obey |ao| > |ai| > |a2|.... If we are only 
interested in the scattering lengths Qq, it is enough to 
calculate the scattering amplitude in the static limit. 



i.e. when s 



4m;, t 



and u 



FIG. 1: Tree level diagrams 



Note that, the linear sigma model is in a better 
agreement with the experimental results than first order 
chiral perturbation theory. 



«o = ^T\s ^ 4ml, t^O,u^O) (9) 



III. PION-PION SCATTERING AMPLITUDES 

The diagrams shown in figure 1, where the continuum 
line denotes a pion, and the dashed line a sigma meson, 
contribute to the tt-tt scattering amplitude. The diagram 
with a sigma exchanged meson has to be considered also 
in the crossed t and u channels. 

From these diagrams it is possible to get the results 
shown in table I. The isospin dependent scattering am- 
plitudes at the tree level have the form 



IV. ONE LOOP THERMAL CORRECTIONS 
FOR SCATTERING LENGTHS 

For our calculation of the thermal corrections to the 
scattering lengths, we will use the real time formalism. 
At the one loop level it is enough to use the Dolan- Jackiw 
propagators |12l |. A general review about the real time 
formalism, beyond the one loop level can be found in [l^ , 
[l^ . In our case the most relevant thermal contributions 
will be related to the pion sector, due to the Boltzmann 
suppression in the case of the sigma meson and/or nucle- 
ons. Therefore, for the pion propagators we will use 



3 



A(A;o, k,m^)Saf3 = ( 



+ 27mB(|fco|) 

+ le 

Sik^-P -ml))6^p, (13) 



-A(fco, fc, m^r) A(fco — 2mT, fc, 771^) (15) 



is the Bose-Einstein distribution. 



where nB{x) 

For the a meson propagator, we wiU restrict ourselves to 
the r = case 



iMfc =16A*i;^( 



4m?. 



d^k 

(2^ 



^ [A(fco, k,ma 



A(fco + m^r, A:, m^) A(A;o ~ "^7r, fc, "t-tt)] 



(16) 



A{ko,k,m^) 



fr2 



(14) 



There are many diagrams that contribute to the pion- 
pion scattering amphtude at the one loop level. For each 
one of these diagrams we have to add also the correspond- 
ing crossed t and u channel diagrams. In figure 2 we have 
shown only the s channel contribution. 

We will give the analytic expression only for the first 
two diagrams (a,b) shown in figure 2. It should be noticed 
that, when it corresponds, symmetry factors, isospin in- 
dex contractions and multiplicity factors should be in- 
cluded. 



In some of these diagrams we have to deal with terms 
S^{k), and In those cases the following identity is 

useful [H 



J-(^)"A 



A 



Tl + l 



(17) 



Actually this expression can be proved in the full 
matrix formalism of Thermofield Dynamics. However, 
for some diagrams, we checked our results using also the 
Matsubara (imaginary time) formalism. 

When dealing with integrals of the following shape. 





I^2mJ dkoF{ko){-^[^^^^^mko - Wk) (18) 

+ 6iko + Wk)) + C-^^^)^mo - wk) + 6iko + wk)]) 

where F{ko) is an arbitrary function and Wk = 
fc2 + , we can integrate by parts getting / = 7^^. + 



, where 



27r hm ( 



^(^)"B(|fco|) 



Awl 



,F{ko)sgn{ko){csch{^-^))^ F{wk)nB{wk) 



16w|T 



(19) 



and 



= 27r lim (- 



»(^W(|fco|) 
Awl 



F{ko)sgn{ko){cschC-^)f F{~Wk)nB{wk) 



l&wlT 



Awl 



)• 



(20) 



FIG. 2: Relevant one loop diagrams 



After taking all the diagrams into account, we find the 
following expression for the one loop thermal corrections 
to the TT-TT scattering amplitudes 
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A{s, t, u)t = -A^(14 + 80A\>^( 



4to^ — ^ 



1 



4m^ — 



+ A^(16w2 + 32A^w''( 



1 

4m^ — ml 



)) 



+ 32 A'^t;^ [E'o (tott , Too- , Wfc , T) + £'o ( - m^r , Too- , Wfc , T)] 
+ 32 A^w'' [i^i (to^ , Too- , Wfc , T) + £'j ( - , mo- , Wfc , T)] 

+ 8AV( ^^^^_^J X [Foi(m.,m.,z.,,r) 



1 



+ ( -TO, , , Wfc , T)] + A^ [96 A^ ^4 ( ) 

4771^ — TOJ. 

..2l 



+ 16v'^]Fl{m^,ma,Wk,T), 



(21) 



s, u)t = -4A*A(wfe, T) - A4(11 + 80X''v^{ j) 

-mj 

+ 48AV(^)2)i?(z«fe, T) + 4A%2(^) 

X ( / , +^)C(u;fe,r) 
TO^ — mj. 

+ 16A''?;^[D^(TO,,mCT,'u;fc,T) + (-771,^, to^t, w^, T)] 
+ l6X^v'^[D'^{m-^,m^,Wk,T) + D^{-mTr,ma,Wk,T)] 

+ A^(32w2+32A\>4(^^)) 

X [EQ{m^,ma,Wk,T) + El{-mT,,m„,Wk,T)] 

+ 16A^w''[i?Q(m7r,TOcr,7«fc,T) + El{-mT,,ma,Wk.,T)] 

+ A^ (48 A2z;4 ( ^ ) + 87-2 ) („,^ , , 7i;fc , T) 

+ Fo^(-TO,,77lCT,Wfc,r)] 

= ^(7i,t,s)T, (22) 

where we have introduced the following definitions 



B{wk,T) 



(27r 
-( 



'-ld\kr-^, (23) 



(27r)3^ 277;^. u., 
d^k UBiwk) 



(2^)3 Wk 

(27r)3 Wk 



), (24) 
(25) 



X '--^(TzrA-r)] (26) 



/"(cc, ?/) 4a;2 + 4x7(;fc 



d^k {x + Wk)nBiwk) 



{27r)^'2wlr+^{x,y)f"^{-x,y) 



^ {csch{^)f ^ riBiwk) ,, 



16r77;2 477;3 ^V"(x,2/)/™(-x,7/) 



)] (27) 



F;^{x,y,Wk,T) 



and where 



d3fc UBjWk) , 1 , 

(27r)3 Wk ^r'{x,y)r^{~x,yy' 

(28) 



r{x,y) = {2x^ + 2xwk-y^y 



V. NUMERICAL RESULTS 



(29) 



The thermal corrections must be calculated numeri- 
cally. The different parameters in our expressions are 
renormalized at T = , since thermal corrections do 
not add new ultraviolet divergencies. The linear sigma 
model, excluding the nuclcons, has three parameters: 
TO^, /tt and A^. The first two parameters, to^ and /tt, are 
given by experiments and the third one is a free parame- 
ter. Notice that f-^ is related to the vacuum expectation 
value V. In fact, at the tree level /,r = v. The three 
parameters are not independent. If instead of f^^ we use 
the vacuum expectation value v and consider a mass of 
the sigma meson ttt-c, — 700 MeV, we have A^ = 7, ti = 90 
MeV; if A^ = 5.6, v = 120 MeV 0]. We know, however, 
that the mass of the sigma meson is about m^ = 550 
MeV [l7[. Therefore, we need to find new values for A 
and V associated to the new lower mass of the sigma me- 
son. We found A^ = 4.26 and 7; = 89 MeV, following 
the philosophy presented in [l^. The authors used the 
Pade approximant method, but they also suggested an 
analytic approach in order to find the variation of one 
parameter if the other two change as a consequence of 
one loop corrections. The values given above for A" and 
V were found following this procedure. 
The scattering lengths, including our thermal corrections 
are given by 



aO(T) = Q^^,T^^Ais,t,^)T + 2Ait,s,u)r ^^^^ 



al{T) = -0.041 + 



327r 
2A(t,s,u)T 
32^? ■ 



(31) 



The behavior of the normalized scattering lengths 
°° and ""''J'' are shown in figure 3. Notice that ah(T) 
vanishes identically. We normalized with the zero tem- 
perature values of the scattering lengths, according to a 
two loop perturbation calculation |18| . It is interesting 
to remark that ag(T) grows with temperature and the 
same occurs with the a§(r). Similar calculations have 
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FIG. 3; Scattering lengths normalized to T = 



been done in the context of the Nambu-Jona-Lasinio 
model [l^ and in the frame of chiral perturbation theory 
[20| . The results from the Nambu-Jona-Lasinio analysis 
do not agree with our conclusions. In this approach 
both scattering lengths remain almost constant, but 
with the same growing tendency. The results from chiral 
perturbation theory agree quite well with our conclusions 
for aQ(r). Nevertheless, a§(r) is almost constant in this 
approach. In \2^, a different determination of ao{T) 
based on the Heat Kernel expansion technique applied 
to the linear sigma model, is presented. We agree with 
their results. 
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